Elementary Thermodynamics of Trapped Particles 
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I develop simple thermodynamic relations for a collection of noninteracting classical particles 
confined in a harmonic trap. The volume of such a trap is not a good thermodynamic variable, so 
conventional expressions of the first law of thermodynamics and the ideal gas law must be modified. 
I use the frequency of oscillations about the minimum of the trap as an external parameter char- 
acterizing the confinement, and derive elementary relations between particle number, temperature, 
energy, oscillation frequency, and a generalized pressure, that are analogous to conventional ther- 
modynamic relations for an ideal gas in a rigid volume. I also discuss heat capacities for trapped 
particles. 

PACS numbers: 05.70.-a, 05.50.Ce, 51.30.+i 



I. INTRODUCTION 

The derivation of the ideal gas law is covered in es- 
sentially every modern text on thermodynamics and sta- 
tistical mechanics. This law applies to non-interacting 
classical particles in a rigid container, and expresses the 
familiar relationship 



PV = NkT 



(1) 



between the pressure P, volume V, number of particles 
N, and the temperature T; k is the Boltzmann constant. 
It is also shown that the energy of of an ideal gas is given 
by 



E 



-NkT. 



(2) 



Recent experiments on cooled neutral atoms (which have 
ultimately led to observations of Bose-Einstein conden- 
sation) have been performed in atom traps which in 
the ideal case have confining potentials of infinite range. 
Thus volume is not an appropriate thermodynamic vari- 
able. It is instructive to consider how the familiar ele- 
mentary relations of Eqs. (|l|) and (||) must be modified 
in these circumstances. 

In this paper I consider atoms in isotropic harmonic 
confining potentials. For particles of mass m and en- 
ergy E confined by a one-dimensional harmonic poten- 
tial with angular oscillation frequency to, the amplitude 
of oscillation is yj2E/(muj 2 ). This suggests using the fre- 
quency of oscillation about the minimum of the trap as 
the externally determined parameter characterizing the 
confinement of the particles in the trap, and in this note 
I consider the development of thermodynamic relations 
which involve 10 rather than the volume V. 

I follow the approach used in several modern texts to 
derive the ideal gas law (see, for example, [Q, |[ [| [I|). 
The canonical ensemble partition function is derived in a 
semi-classical manner, i.e., using information about the 



spacing and degeneracy of quantized energy levels, and 
the partition function is then used in combination with 
the first law of thermodynamics to derive relationships 
between the thermodynamic variables. 



II. FUNDAMENTAL EQUATIONS 

The first law of thermodynamics is an articulation of 
the work-energy theorem. For a gas in a volume V this 
is expressed as 



AE = AQ - PAV. 



(3) 



Although the volume isn't a relevant parameter for par- 
ticles in a trap, the confining potential can do work on 
the particles when it is altered. We can express this idea 
in a modification of the first the first law, 



AE = AQ + VAuj, 



(4) 



in which V is a "pressure" for which I will derive an 
expression below. (I use the term "pressure" loosely; this 
quantity does not have the dimensions of force per unit 
area.) Note that I have chosen the sign of the second 
term on the right side of Eq. (^ to be positive, reflecting 
the fact that an increase in u> corresponds to an increase 
in the strength of the confinement, which increases the 
density of the particles and effectively compresses them. 

The conventional arguments leading to the definition 
of the Helmholtz free energy 



F = E -TS 



(5) 



are not affected by the change from confinement in a 
rigid volume to confinement in a trap. Combining the 
modified first law with the definition of free energy leads 
to the relationship 



dF = -SdT + Vduj, 



which implies 
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V = 



dF 
du> 



(6) 



(7) 



T,N 



2 



and 



(8) 



The relationship between the Helmholtz free energy 
and the partition function is based on general arguments 
regarding entropy, so that for particles in a trap it is still 
true that 



F(T,w,N) = -kT\nZ(T,uj,N). 



(9) 



III. 



EXPLICIT PARTITION FUNCTION AND 
CONSEQUENCES 



The energy of the particles is 

'dlnZ x 

W / 



E = - 



= 3NkT, 

which combined with Eqs. (|l3|) shows that 
E = Puj. 



(15) 



(16) 



The linear relationship between energy and oscillation 
frequency is to be expected given the fact that energy of 
each quantized single-particle energy level is of the form 
mhco (where m is an integer). 



The partition function for a single particle in a a har- 
monic trap is 



Zi = 



-0ei 



f(e)e-^de 



2(hujf 



kT 



(10) 



In the derivation above I have assumed that kT ^> %lo 
and converted the discrete sum into an integral, and I 
have used the density of states /(e) appropriate for a har- 
monic oscillator potential. Using standard arguments, 
the partition function for N non-interacting particles in 
a dilute gas is 



(s) 



N (kT 



(11) 



where in the last line I have used Stirling's approximation 
for large N. 

The Helmholtz free energy is thus 



F 



-kTlnZ 



= -NkTln 



kT 

N V hw 



(12) 



The "pressure" P is given by 
V = 



dF 
du 
3NkT 



T,N 



(13) 



which gives the analog to the ideal gas law for the trapped 
particles: 



Plu = 3NkT. 



(14) 



IV. COMMENT ON HEAT CAPACITY 

For conventional ideal gases the heat capacity at con- 
stant pressure is 



(17) 



and this is related to the heat capacity at constant vol- 
ume Cv by the familiar equation 



C P -C V =P 



dV 
df 



Nk. 



(18) 



For trapped particles the heat capacity at constant lo 



is 



O = [ ?jf) ■! V/,-. 



(19) 



This makes sense because the translational motion of 
each particle contributes 3fc/2 to the heat capacity, and 
the potential energy of the three-dimensional harmonic 
oscillator confining potential contributes an additional 
3k/ 2. The difference between the heat capacities at con- 
stant lo and at constant V is 



= -P(^ 
\dT 

= -3Nk. 



This means that 



C v = 0. 



(20) 



(21) 



In other words, the energy needed to raise the tempera- 
ture of the particles by AT all comes from the work done 
on the particles by the increase in to that is necessary to 
keep V = E/uj constant. 



V. CONCLUSION 

For particles trapped in smoothly varying long-range 
potentials, volume is not an appropriate thermodynamic 
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variable. Thermodynamics can, however, be developed 
using other parameters characterizing the confinement of 
the particles. In this paper I have considered particles in 
isotropic harmonic traps, and I have used the frequency 
of oscillations about the minimum of the trap as the con- 
finement parameter analogous to V. I have derived the 
analog to the ideal gas law for this simple case, and also 



formulas for appropriate heat capacities. These simple 
relations can be used to solve many thermodynamic prob- 
lems for trapped particles that are analogous to problems 
for gases in containers with rigid walls that are posed in 
introductory texts. Generalizations to anisotropic traps 
and more complicated potentials are certainly possible, 
and would make good student projects. 
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